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A Contribution to the Thermodynamical Theory of Explosions. — 

Parts I and II. 

By Sir J. B. Henderson, D.Se. and Prof. H. E. Hasse, D.Sc. 
(Communicated by Sir J. A. Ewing, F.R.S. Eeceived May 6, 1921.) 

Introduction. 

The following paper deals with the application of the fundamental ideas of 
modern physical chemistry to the science of internal ballistics, and aims at 
eliminating some of the empirical methods and the arbitrary assumptions 
which have been necessary in the practice and theory of that science in the 
past. The object of the science is the calculation of the pressure and 
temperature of the products of explosion throughout the motion of the 
projectile along the bore of the gun, both during the burning of the propellant 
and after the burning is completed. It is necessarily a branch of physical 
chemistry, since it implies a knowledge of the composition of the products of 
explosion and the variation of these products with temperature, also a know- 
ledge of the variation of the specific heats of gases with temperature. Hence 
every advance in our knowledge of these branches of physical chemistry 
will enable a further step to be made in the science of internal ballistics. 

To cover our ignorance in this respect it has been customary in internal 
ballistics to make arbitrary assumptions as to the connection between pressure 
and volume behind the projectile in the gun. The assumption made for 
example by Mansell and Hadcock* is that the equation, p(v—h) n /z n = constant, 
is satisfied, the value of n being different during burning and after burning is 
completed, z being the fraction of the charge already burnt. 

The present paper makes no assumption of this kind, but contains an 
application to the problems of internal ballistics of the recent advances in the 
thermodynamics of dissociation and in the thermodynamics of the variation 
of specific heats of gases with temperature, with the result that the whole 
indicator diagram of the gun can now be calculated from the chemical 
composition of the propellant, combined with a table of heats of formation 
from the elements of the propellant and its products, and one set of labora- 
tory experiments to determine the rate of burning of the solid propellant 
and its variation with pressure. Physical chemistry is not yet sufficiently 
advanced to enable this rate to be calculated. 

Part I contains the calculations of the maximum temperature and pressure 

*■ A good review of previous work is given by Sir George Hadcock, F.E.S., in his 
paper on Internal Ballistics, 'Roy. Soc. Proc.,' A, vol. 94 (1917). 
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for an explosion in a closed vessel, also the calculation of the curve of 
reversible adiabatic expansion of the products and incidentally the ideal 
indicator diagram. This part was written in 1913-14 and publication was 
then prohibited. 

Part II was written in 1915-16 and contains the calculations for the gun, in 
which the variation with time of the volume of the chamber, of the com- 
position of the products and of the temperature and pressure of explosion, are 
all worked out point by point and the value of the maximum pressure is 
determined. 

This complete theory was put to very severe test by application to many 
practical problems during the War, made by that distinguished scientific 
artillerist, Captain Maxwell Lef roy, R.N., member of the Ordnance Committee, 
with very satisfactory agreement between theory and practice. Incidentally 
it may be noted that since actual direct experiments on specific heats of gases 
have oeen limited to temperatures below 1500° C. and pressures of a few 
hundred pounds per square inch, whereas in ballistics temperatures of 3500° C. 
and pressures of 20 Tons/in. 2 are common, the application to ballistics of the 
theoretical thermodynaniical equations forms a severe test of the theory. 

Part I. — The Calculation of the Temperature of Explosion, the Curve of 
Acliabatic Expansion, ivith Special Reference to the Explosion of Cordite and 
the Ideal Indicator Diagram* 

By Sir J. B. Henderson, D.Se. 

The general theory will be illustrated by an example of its application to 
the explosion of Mark I. cordite, and the results for M..D. cordite are added, 
but it will be understood that the method is equally applicable to any 
propellant. 



1 grin, of 
Mark I. 
cordite 



"37 grm. guneotton 



n.Kc te T l • J 0*0800 „ hvdrogen . 

0*58 „ nitroglycerine = < n . 1eKO ■•*. to =s 

~ rtr " ,P J I •1553 „ nitrogen 

'Ob ,, vaseline L irQ9 h to 

" l_0 5837 ,, oxygen 



"0 '2308 grm. carbon fO -01924 grin, atom of C. 

0*0300 „ „ H. 
O'Olll „ „ N. 
0*0365 „ „ O. 



The products of explosion will contain C0 2 , H 3 ,- CO, H 2 Q and N 2 with 
possibly small quantities of NH 3 , CH 4 , ISTO ; it is also possible that some of the 
CO2 and H 2 are dissociated, but whatever the composition, it will be assumed 
that at any temperature the composition represents a state of thermodynamic 
equilibrium so that the equations for such a state are applicable. We shall 
first of all neglect the small quantities of NH3, CH 4 and NO present, also the 
dissociation of CO2 and H 2 0, and prove later that these are negligible. The 

* The author is much indebted for valuable assistance and advice on the chemical side 
to Dr. Thomas Ewan, of Glasgow. 
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gaseous reaction which varies with the temperature is then the water-gas 
reaction 

C0 2 + H 2 :Z^CO-hH 2 0. 

Let p, q, r and s be respectively the numbers of gramme-molecules per 
gramme of products, and let C, 0, H and N" be the numbers of gramme-atoms 
of C, 0, H and 1ST respectively in the products, then 

v + q = G; 2p + q--=: 0; 2r+2s = H; ^ = K* 

M J " ' pv 

where K = '- G0 J x [ H »°3 

[CO a ] X [H 3 ] ' 

Values of log K c for the water-gas reaction may be obtained from works on 
physical chemistry. 

Since K c varies with temperature, it is necessary to assume first of all a 
value of K c . then solve for p 9 q y r and s, and if the values so obtained lead to a 
value of the temperature considerably different from that associated with the 
assumed value of K C} then the calculations have to be repeated with a new 
corrected value of K c . A large error in K c , however, makes a small difference 
in the products, 

Taking K c = 14 which corresponds to a temperature of 3400° abs. 

p = 0*00486 grm.-mols. per grm. cordite. 
q == 0*01438 
r = 0*00262 



>) y> 



» ?; » 

s = 0-01238 



N 2 = 0-00554 



'J 5> !•) 



>' J' 5) 



Sum = 0*03978 

The characteristic equation for a gas, or mixture of gases, at high tem- 
peratures andjpressures \%p{y~ b) = Rd, where v is the specific volume of the 
products, h the specific volume of the molecules themselves, which may 
be taken as the specific volume of the solid explosive. The value of E 
for a mixture of gases is the sum of the concentrations of the constituent 
gases in gramme-molecules per gramme, multiplied by a constant which is 
2760 ft Lbs./ lb. x deg. C. Hence 

E = 0-0398 x 2760 = 109*8 -J±- Lbs ' = 3*25 atmos ' x cm ' 3 



lb, x deg. O. grm. x deg. 0.* 

Thermal Values. 

Berthelot has given the heats of formation from the elements of collodion, 
guncotton, nitroglycerine, vaseline and all the above gases formed in the 
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explosion, all at 15° C. The heat evolved per gramme of cordite at 15° C. is 
next calculated from Berthelot's figures combined with the above values of 
p, q, r, etc., which corresponds to the heat evolved in the explosion, or the 
heat evolved if the gases were cooled from explosion temperature to 15° G. 
without altering their composition. 



Molecular 
\ weiglifc. 


Heat of formation 

per grm.-mol. 

from elements 

(at 15° 0.). 


Heat of formation 
per grin. 


Collodion, Co^isiOsgNy 

Gruncotfcon, C24II09O42N] t 


1053 
1143 

227 

282 

44 

28 

18 

18 

* 


cals. 
696,000 
624,000 
98,000 
97,100 
94,300 
26,300 
58,000 1 
68,400 J 


660 
546 
432 
344 


Nitroglycerine, O3H5N8O9 

Vaseline, C20H49 (?) 


Carbon dioxide, CO.-. 


7 4hf 

Carbon monoxide, CO 

Steam, H<>0 





Since the heats of formation are determined at atmospheric pressure, it is 
necessary, in order to get the heat of formation Qi at constant volume, to add 
the energy which has been expended in doing external work, namely, 
550 calories per gramme-molecule of gas produced. Instead of adding this 
to each item, it will be added at the end of the heat balance. 

Hence the heat balance of the explosion of cordite, Mark I, stands as 
follows, assuming guncotton to be a mixture of 70 per cent, guncotton, 
C24H29N11O42 and 30 per cent, collodion C24H31N9O38 :— 

Heat absorbed by the decomposition at 15° C. of — 

07 x 0*37 grm. guncotton = 546 x 0*259 = 141*4 cals 
0*3x0*37 „ collodion =660x0*111= 73*3 
0*58 grm. nitroglycerine = 432 x 0*58 = 250*5 
0*05 „ vaseline = 344 x 0*05 = 17*2 



>? 



?> 



Sum = 482*4 cals. 



Heat evolved by the formation at 15° O. of— 

0*00486 grm.-mol. OO a = 0*00486 x 94300 

0*01438 „ CO = 0*01438 x 26300 

0*01238 „ steam =0-01238x58000 
0*0398 



>) 



of gas = 0*0398 x 550 



Sum 
Hence the total heat evolved per gramme of cordite 



458'4 cals. 


378-2 „ 


718-5 „ 


Zt x. \j ,, 


1577*0 cals. 


1094*6 cals. 



On the Thermodynamical Theory of Explosions. 465 

The heat evolved in cooling the gases from explosion temperature to 
15° C. must be equal to 1094*6 calories, whence the temperature of explosion 
is obtainable if the specific heats of the products are known. 



Specific Heeds. 

Experimental determination of the specific heats of gases by direct 
measurement is limited to temperatures below 1500° C, and without a 
thermodynamic theory of the variation of the specific heats with tem- 
perature, further advance in the theory of explosion would be difficult, if 
not impossible. Such theories are those proposed by Einstein and by Nernst 
and Lindemann, based upon Planck's quantum theory. 




WOO 1200 1400 1600 /800 2O0O 2200 2400 2600 2800 JOOO 

Fin. 1 



The curves in fig. 1 show the variation of the specific heats of C0 2 , H 2 0, 
N, and H, with temperature as given by Bjerrum, calculated by the 
Nernst-Lindemann formula 



dE E f" x 2 e 



X 



-- < 



--f- 



(x/2Y 



2 e x/2- 



d0 2 L(^~-l) 2 (e*l 2 -Vf 



where x 



hn 

w 



The following table is prepared from these curves :- 



vol. c. — A. 



2 H 
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Grm.-mol. of products 

per grm. of 

cordite. 


t = 2000° 0. 

Mol. i MoLhoafc 
heat. . ^m.-mol. 
i oi products. 


Mol. 

heat. 

10 -57 
5*94 
5*78 
9-46 


2500° 0. 

Mol. heat 

x grm.-mol. 

of products. 

'0514 
-1184 
'0151 

0*1171 


t = 

Mol. 
heat. 


3000° C. 

Mol. heat 

x grm.-mol. 

of products. 


C0 2 = 0-00486 

00 + 1*3 - 0-01992 ... 

H 2 - 0-00262..... 

H 2 -0-01238 ......... 


10 ' Oa 

5-77 
5'62 
8-42 


'0501 
1150 

-0147 
0-1042 


10*75 
6-15 
5-95 

1.1 -02 


'0523 
0-1226 
-0156 
0*1365 


Heat capacity of products — 


-2840 




'3020 


— _ 


-3270 



Now, the product heat capacity x temperature = heat contents, and the 
temperature has to be calculated at which this product is equal to the heat 
evolved in the explosion, namely, 1094*6 calories per gramme of cordite. 



| Temperature. 


Heat capacity. 


2000 
2500 
3000 


*2840 
-3020 
'3270 



Product. 


Difference. 


568 
755 
981 


187 
226 



Plotting the product on a temperature base and reading off the tem- 
perature corresponding to a product 1095, we find 3210° C, or 3483° 
absolute. 

In Service Chemistry a simple approximate formula is given for Mark I. 
cordite : — C5H8N3O10. If the above claculations are repeated, using this 
formula the results are :—p = 0*00572, q = 0*01278, r = 0'00202, s = 0*01278, 
N 2 = 0*00556 and E = 3*175. The heat evolved by the formation of the 
products = 1638-2, and, assuming the heat absorbed by the decomposition of 
the cordite molecule the same as above, viz., 482, the heat balance is 1156*2 
and the temperature of explosion is 3600° C. absolute, instead of 3500° from 
Berthelot's composition formulae. 

The calculations have also been repeated for MJD. cordite, containing 
65 per cent, guncotton, 30 per cent, nitro-glycerine, and 5 per cent, vaseline, 
with the following results : — p = 0*0031, q = 0*0184, r = 0*0046, s = 0*0111, 
T$ 2 = 0*0050, Sum = 0*0422, E = 3*44 (atrnos. x cm. 3 )/grm. x deg. C). The heat 
of formation of 1 grm. of M.D. cordite is 524, and the heat balance works 
out at 913 calories evolved per gramme; and the temperature of explosion is 
2870° C, or 3143° O. absolute, or 340° (I lower than for Mark I. 
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Influence of Initial Temperature of Cordite. 

The influence of the initial temperature may be determined as follows : — 
Let the explosion take place at different initial temperatures, t\ and t 2) the 
corresponding temperatures of explosion being 1\ and T 2 . It is easy to 
imagine the chemical change taking place at ti until the products are the 
same as exist at Ti, and that the products are then heated up from t\ to Ti 
by the heat evolved. Secondly, imagine that the cordite is heated up from 
t\ to t 2 by adding heat h from external sources, then let the chemical change 
take place at t 2 until the products are those existing at T 2 , the heat evolved 
being then used to raise the temperature of the "products from t 2 to T 2 . If 
the heat h were not supplied from external sources, the final temperature in 
the second case must be Ti, since the temperature of explosion is independent 
of the path followed in the explosion. Hence the heat h is equal to that 
required to raise the temperature of the products from Ti to T 2) neglecting 
the small difference in the products at these temperatures, so that 

(specific heat of solid cordite) {fa—ti) 

= (specific heat of products at constant volume at 
temperature of explosion) (T 2 — Ti). 

Now, the specific heats of solid cordite at ordinary temperatures are 0*365 
for Mark I. and 0*347 for M.D., as given by Sir Eobert Eobertson, and, for 
the former, the specific heat of the products at the temperature of explosion 
is 0*372, as calculated from the above theory. Hence the rise in the 
temperature of explosion is 0*365/0*372 of the rise in the temperature of 
the solid cordite, so that the temperature of explosion is practically increased 
by an amount equal to the rise of temperature of the solid, cordite. For an 
increase from 60° F. to 80° F., the temperature of explosion is therefore 
increased by 20° F., or 11° C, which is within the limits of error in the 
determination of the temperature of explosion itself. 

Pressure of Explosion. 

The temperature of explosion having been determined, the pressure for 
any particular density of loading follows at once from the characteristic 
equation. For Mark I. cordite, o = 0*641 cm. 3 /grm., and, if the density of 
loading be 0*5 grm. per cm. 3 of bomb capacity, 

p = 8370 atmos. = 55*1 tons/in 2 . 



2 h 2 
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Correction of Results for Dissociation. — Dissociation of CO2 at 3500° abs. and 

Pressure of Explosion. 

logwK, = log 10 [C S 3 ' [ y = -1-6 = 2-40. 

[002/ 

With a density of loading 0'5 grm./em 3 , or 1 grm. per r 002 litre, the con- 
centrations are at once determined by the composition of the products on 
p. 463. 

[C0 2 ] = °^^ = 2-43 grm.-mok/litre, 
rpm _ 0-01438 _ 7 

[G0] - IFooT ~ 7 19 



J5 )> 



whence [0 2 ] = 0*00287, or 0*0012 grm.-mols. of 2 per gramme-molecule of 
C0 2 , and since the dissociation of 2 grm.-mols. of C0 2 give 1 grm.-mol. of 2 , 
the percentage of C0 2 dissociated is 0*24. If the density of loading were 0*05 
the percentage dissociated would be 2*4. 

Dissociation of Steam. — Dealing with the H 2 in a similar manner and 
taking log K c = 5*97, gives 0*06 as the percentage H 2 dissociated at a density 
of loading of 0*5, which is negligible. 

In a similar manner we get [CH 4 ] = 1*3 x 10~ 6 which is negligible. CH 4 is 
formed in cooling the products. 

[NO] = 0*0160 grm.-mols. per litre. 

= 32xl0~ 6 grm.-mols. per gramme cordite. 
= 0*96 per cent, of the products of decomposition. 
[NH3] == 0*41 grm.-mols./litre. 

= 0*00082 grm.-mols. per gramme of cordite. 
= 1*4 per cent, of the products of decomposition. 

All these are very small as compared with the main products, C0 2 , CO, H 2 , 
H 2 0, and N2, and are negligible where an accuracy of about 1 per cent, is aimed 
at. If any of them were not negligible, it would be easy to allow for it by. 
successive approximation. 

Heat Contents of the Products of Explosion and Curve of Adiabatic Explosion. 

The composition of the products, if in chemical equilibrium, depends only 
on the temperature, and the intrinsic energy in each of the constituents 
depends also only on the temperature, hence the heat contents of the products 
of explosion is a function of the temperature only. It is therefore necessary 
to take an arbitrary datum both for temperature and for composition of 
products. The datum taken for temperature is 0° C. and the datum for 
composition is conveniently taken as the composition at 1500° C. abs. 



On the Thermo dynamical Theory of Explosions. 469 

The heat contents is, therefore, the heat evolved in cooling the products to 
1500° abs., allowing combination to take place, and then cooling the products 
further to 0° C. without allowing further combinations to take place. It will 
be more convenient in tabular work to cool the products to 0° C. without 
allowing combination to take place and then allow the product to combine at 
0° O. until the composition is that corresponding to 1500° abs. The result 
must be independent of the path followed, since it is a function of the 
temperature only. 

In Table I the composition of the products of Mark I. cordite is given at 
1500°, 2000°, 2500°, 3000°, and 3500°. abs., each column being calculated by 
the method given on p. 463. The heat of reaction for water-gas is added, and 
these figures might be used if the products were imagined to* combine at 
these temperatures to the standard composition before being cooled to 0° C. 
They are added chiefly to show that a maximum heat of reaction exists at 
about 2000° abs. This maximum must have a considerable effect in influencing 
the final composition of the products formed in bomb calorimeters. The heat 
of reaction at 0° C. is 10,420 calories per gramme-molecule, and its value at 
any other temperature t is obtained from the equation 

Q + 1 (mean specific heat of products) 

■ = Qo + ^ (mean specific heat of the combining gases). 



Table I. — Composition of the Products of Combustion of Mark I. Cordite at 

various Temperatures. 



Temp. deg. G 


1227° 
1500° 


1727° 

2000° 


' 2227° 
2500° 


• 

2727° 
3000° 


3227° 
3500° 


Temp, abs 




K c 


3-71 


6 -61 


9-5 


12-6 


16-2 




0O 2 


-00700 
-01224 
-00476 
-01024 
-00554 


0-00600 
-01324 
O -00376 • 
-01124 
-00554 


-00542 
-01382 
-00318 
0-01182 
-00554 


-00502 
O -01422 
-00278 
-01222 
-00554 


-00467 
-01457 
• 00243 
01257 
-00554 


GO 


H. 


it 

Steam 


N 2 






O -03978 


-03978 


-03978 


-03978 


O -03978 




Heat of reaction 


13,342 


14,082 


13,717 


11,866 


6,548 



Table II contains the results of the calculation of the heat capacity of the 
products of explosion of Mark I. cordite, the datum of reference being the 
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products at 1500° abs., cooled to 0° C. The total heat contents given in the 
last line of the Table are represented within 1 per cent, by the equation 

H = 15*2 + 01223 + 55+ xlO~ 6 3 

where H is in calories per gramme of products of explosion. 
Table III corresponds to Table II, but refers to M.D. cordite, 

A diabatic Expansion, 

The equation for a reversible adiabatic expansion expresses the fact that 
the loss of intrinsic energy is equal to the external work done, or 

p&v= -(||W= -(0-1223 + 110-8 xlO~~ 6 0)S0, 



and miGep(v—b) = K0 it follows 

civ cW 



v-b ~ 110 



(0-1223 + 11 0-8 x 10 ~«0) 



, 1 v--b 0-1223 ! 0i , 110-8 x 10~ 6 /V) m 

whence log, == — - — log tf — + — (0i — 0) } 

v\ — b li u li 

and substituting the values for K we get 

logio— j = 1'552 log 10 ^ + 617x 1O-«(0i-0), 

which is the equation for adiabatic expansion. It is easily shown that this 
equation cannot be represented with sufficient accuracy by an equation of 
the form p(v — b) n = constant, which necessarily involves H being a linear 
function of 0, whereas it is a quadratic function. This latter type of equation 
could be used for small changes of temperature, n varying with the tempera- 
ture according to 

II dll = 0-1223 + 110-8 xlO- 6 0. 



n — 1 d0 

Remits for M.D, Cordite. 
H = 14-68 + O\L3230 + 52-57xlO- 6 3 . 
E = 3*486 x 10~ 6 ergs./grm. x deg. 0. 
= 0*0831 grm. cals./grm. x deg. 0. 
Adiabatic expansion is given by 

l ogl0 2=L = 1-593 log § + 549 x 10~« (0i-0). 

V\ — u 

Heat balance = 913 cals./grm. ; temperature of explosion = 2870° 0. or 
340° lower than Mark I. 
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Ideal Indicator Diagram and Mechanical Efficiency of a Gun. 

The ideal indicator diagram may be considered as that obtained by a sudden 
explosion of all the propellant before the projectile begins to move, to be 
followed by an ad.iabatic expansion of the products. The maximum pressure 
obtained in this way is much greater than the gun could stand ; hence it is 
fairer to take the ideal as that obtained by burning the propellant at the 
maximum working pressure for which the gun is designed, to be followed by 
adiabatic expansion. 

Fig. 2 shows both types of ideal indicator diagram for MJD. cordite, and, in 
addition, the curve of pressures for irreversible adiabatic expansion, as in 
throttling, in which no work is done, corresponding to explosions in closed 
vessels of different capacities ; also the ideal indicator diagram assuming an 
elliptic curve of rise of pressure; also the uniform pressure at which the 
cordite would just be all burnt at the muzzle of a gun having a ratio of 
expansion v/v equal to 4*09. 

Mechanical Efficiency of the Gun. 

Taking the mechanical efficiency of a gun as the ratio of the sum of 
the kinetic energies of the projectile as it leaves the muzzle, due to transla- 
tion and rotation, to the work represented by the area of the ideal indicator 
diagram, as calculated above, the following Table is evolved: — 



Gun. 



Table of Efficiencies (M.D. Cordite). 



Density | Mass 

of | of 
loading, charge. 



6" Mark VII ... 

7 -5" Mark III... 

9 *2" Mark XI 

12" Mark X 

13 *5" Mark V ... 



0-463 
0*404 
0-396 
0-401 
0-419 



lbs. 
28f 
54i 

128| 

260 

297 





Ballistics 










(experimental). 




Total 








Thermo- 


energy 




v/v . 






dynamical 


of 


Efficiency 




Maxi- 


Muzzle 


work. 


pro- 






mum 


Telocity. 




jectile. 






pressure. 


V 










Tons/ins. 2 


ft./sec. 


ft. Tons. 


ft. Tons. 


per cent. 


4-986 


16-55 


2761 


6,237 


5,315 


85-2 


4 -993 


17-79 


2781 


12,400 


10,790 


87*0 


3 -990 


17*76 


2870 


25,990 


21,820 


84-0 


3 954 


16*29 


2735 


50,630 


44,340 


87*6 


4-850 


18-78 


2501 


67,610 


60,970 


90*2 
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RATIO OF EXPANSION 
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PART II. — The Bate of Bise of Pressure of the Explosion of Cordite in a Gun. — 

The complete Indicator Diagram of the Explosion. 

By Sir J. B. Henderson, D.Sc, and Prof. H. E. Hasse, D.Sc. 

The object of this part of the paper is the application of the results of the 
thermodynamical theory of Part I to the general problem of Internal 
Ballistics, the variation of- temperature and pressure behind the projectile as 
it moves along the barrel of a gun. The mathematical theory which follows 
differs from that given in work previously published,* in that it is more 
strictly in accordance with thermodynamical theory ; Charbonnier's equations 
form a special case of the theory, but are not treated in the same way as in 
the present paper. 

The essential features of the theory of Part I of this paper are (1) the 
determination of the temperature 6 of the explosion ; and (2) the expression 
of the heat contents, H, of the products of explosion in the form a'-f /3'0 + y'O 2 . 
The loss of heat contents per unit mass of gas, which measures the work 
done per unit mass, and is denoted by h, may therefore be written in the 
form 

h = A + y^-^^T^ 2 . (1) 

For many purposes it is more convenient to -express o -~- 6 as a quadratic 
function of h. It is found from (1) that for M.D. cordite the equation 

<9 o -0 = 3-381 A + 20-45 x W~*h 2 (la) 

will give results in agreement with (1), well within the possible error of (1) 
as deduced from the theoretical considerations of Part I. 

For temperatures differing from that of explosion by less than about 300°, 
h may be expressed as a linear function of 6 in the form 

h=f(d,-d), (iv) 

where for M.D. cordite 

ft. Tons 



/ = 0-28245 



lb. deg. C. 



Fundamental Equations of the Theory. 

In formulating these equations, it should be noted: (1) that all forces 
opposing the motion of the projectile have been neglected; (2) that no 
attempt has been made to allow for the loss of heat to the walls of the gun ; 
and (3) that it has been assumed that the rate of burning of the propellant 
depends only on the pressure of the surrounding gas, and is, moreover, 
proportional to the pressure. The fundamental gas equation is 

p (V—Wi&) = llm0, (2) 

* E.g., the standard work of Charbonnier, ' Balistique Interieure, 3 1908. 
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p and 6 being the pressure and temperature of the gas, which has a mass m } 
V is the volume of the space behind the projectile, m\ is the total mass of 
the charge, and b is the reciprocal of the density of the so Lid explosive. 
The theoretical deduction of this equation depends essentially on this 
interpretation of the constant b, but the equation may be looked upon as 
an empirical generalisation of the similar equation, deduced from experi- 
mental results in closed vessels, and the value of b chosen in accordance 
with these results** 

If M is the mass of the projectile, u its velocity, A the area of the cross- 
section of the bore, including that of the grooves, its equation of motion is 

M f f = P A. (3) 

The scheme of theoretical equations is completed by the equation 

pdV = cl (nth), (4) 

which expresses the work done by the gases on the projectile in terms of the 

total loss, mh, of heat contents at any instant. From these equations, we 

deduce from (2) and (4) that 

-r, dU d(zh) /KX 

U zd 

where z = m/m, h and is therefore the fraction of the charge burnt at any 
instant, and U == Y—mib, which may be called the " free volume " behind the 
projectile. 

Again from (3) and (4) it follows that 

mizh = I Mm 2 } (6) 

the equation of energy, as is otherwise obvious. 

It is known that any propellant when once ignited is burnt in parallel 
layers, so that the rate at which the thickness diminishes is, in accordance 
with the statement made above, given by 

| = /W, TO 

where e is the thickness burnt off from any surface in time t. According to 
Petavelf and Mansell,* /(.p) for cordite is equal to a + aip. Using the 
values given by them, the first term is only of importance for low values of 
the pressure, i.e., in the initial stages of the explosion, or when the densities of 
loading are low, when there are other factors of much more importance to be 
taken into consideration, so that in what follows we shall take f(p)z=ap&s 

* Charbonnier, op. cit., pp. 9-17. 

t Petavel, < Phil. Trans.,' A, vol. 205 (1906). 

J. Manseli, < Phil. Trans.,' A, vol. 207 (1908). 
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being more applicable to normal gunnery practice, leading to a considerable 
simplification in the mathematics of the theory. This is also taken to be the 
law of burning for French propellaifts by Charbonnier. There is, finally, the 
relation connecting the thickness e with the fraction z of the gas burnt off, 
which depends upon the geometry of the propellant, and will be discussed 
later. 

Tt is of interest to consider the connection between this theory and that of 
Charbonnier. Taking (2) and (6), and considering the particular case in 
which the loss of heat contents, h } is proportional to the drop of tem- 
perature, as in (lb), we have, on. elimination of h and 0, 

f 

JMu 2 + *~~ p (V — m\b) — myfOoz 

which is equivalent to Charbonnier's third fundamental equation,* differing 
from it only in the interpretation of the constants involved. As obtained by 
Charbonnier, this equation appears to be an empirical modification of the 
corresponding equation applying to the case when the cordite is. all burnt 
before the projectile moves, and it is curious that it should turn out to be a 
special case of the theory of this paper. 

Derived liquations for Initial Pressure p . 

If it were possible to obtain an expression for the resistance encountered 
by the projectile during the engraving of the driving band, it would not be 
difficult to incorporate it into the general theory of this paper, though the 
resulting differential equations might have to be solved arithmetically. In 
the absence of such an expression, the influence of the driving band may be 
discussed by studying the effect produced on the indicator diagram by raising 
the initial pressure before the projectile begins to move.f 

Suppose therefore that the pressure rises to a value p before the projectile 

moves, a mass m of propellant having been burnt, the equation connecting 

these two quantities being 

p Q = K^o^o/Uo 

as in the explosion in a vessel of constant capacity V , where Uo = Vo—mA 
and V is the volume of the chamber of the gun. Taking f(p) = ap, we have 
from (3) and (7) on integration 

Uu = - («-« ), (8) 

a 

where e is the thickness burnt before the projectile moves. 

* Op. cit., pp. 141, 142. 

t This method has been used by Charbonnier and others. 
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From (6) and (8) we then have 

7 __ ! A 2 ( x 2 . „. i A 2 £i 2 / e e \ 2 
miM& 2 miM<x 2 \#i 01/ 

where 2e\ is the total thickness to be burnt through before the propellant is 
converted into gas. 

Writing B = \A?e^jm{M.a 2 which is the fundamental constant of the theory, 
and putting y = eje h so that y is the fraction of the thickness burnt at any 
time, the last equation reads 

*fc = B(y-y„)» (9) 

The geometry of the propellant gives the relation between % and y in the 
form 

Z=$(y). (10) 

Finally from (5) we get 

jy dTJ __ d (zh) _ 2"B(y—yo) dy ,-..., 

where we have put # = %(h), being the general form of an equation of the 
type (la), and h is further a function of y from (9) and (10). 

The integration of (11) and its reduction to a form suitable for numerical 
calculation is the sole mathematical difficulty, and an accurate analytical 
result can only be obtained in very few cases. For many purposes it will be 
sufficient to give to the temperature a mean value 6, which may be taken 
as the arithmetic means of the temperatures corresponding to the initial and 
final values of y. This gives an approximate integration of (11), which gives 
very accurate results up to the point of maximum pressure, since the variation 
of is small compared with that of y t in the form 

E log (U/Uo) - ~ (Fi-y F a ) (12) 

8 

where ^(^[^^ ™ = i*k 

The form of $(y) has been given for all the ordinary shapes by various 
writers, and is in general a cubic function of y ; (cf., e.g., Ingalls, ' Interior 
Ballistics/ and Charbonnier's work already quoted). In the case of the multi- 
perforated American nitro-cellulose there is a difficulty when the grain breaks 
up into the so-called " slivers," which are cylinders whose cross-sections are 
composed of arcs of circles, and it is not easy to treat these by accurate 
mathematics. Fortunately, in general about 85 per cent, of the charge is 
burnt when this break-up occurs, and approximate methods of calculation can 
be applied to the slivers. In order to avoid any discontinuity in the rate at 
which the pressure changes, the initial volume and surface of the slivers 
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should agree with the initial volume and surface of the more tractable 
geometrical form, by which the slivers are replaced *f or the purpose of the 
calculation. This makes <j>(y) and its differential coefficient continuous at the 
point where the form of the function changes. 

Maximum Pressure of Explosion. 

One of the most important results is the value of the maximum pressure, 
since this can be determined experimentally, and the comparison of the 
theoretical and experimental values is of the greatest importance in testing 
the theory. 

By differentiating (2) and using (2) and (4), it will be found that the 
pressure is a maximum when 

R |w = !<*>• < 13 > 

This result can be simplified if we use (lb) instead of (la) since we expect 
only a small drop of temperature while the charge is burning rapidly before 
the maximum pressure is reached. In this case the left-hand side of (13) is 



k|*(0 o -A//) = R*o-||W, 



so that (13) reduces to 



Ii/0 o = (K+/)Jr(sA), (14) 



or finally, from (9) and (10) to 

K/0 o ^ = 2B(R+/)(y-jr„). (15) 

This equation, in general a quadratic in y, having been solved, the corre- 
sponding value of h is found from (9) and (10), whence 8 can be calculated 
from (16), while the corresponding value of U is found from (12), so that the 
maximum pressure is readily determined. 

So far the investigation applies only up to the instant at which the cordite 
is all burnt. After this there is only the expansion of the whole mass of gas 
to be dealt with and the equation for adiabatic expansion obtained in Part I 
applies ; but if h is used as independent variable, the equation of adiabatic 
expansion is from ^5), since z is now unity, 

p d\J _ dh ___ dh 

the integration of which can be readily adapted to numerical work iu the 
case of cordite, in which is given in terms of h by (la). 
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Results for M.D, Cordite. 

In the case of the standard British propellant, cordite, made up into 
cylinders assumed so long that the burning at the ends can be neglected, 
<t>(y) = 1— ?/, and (15) is therefore a simple equation for the value of y at 
maximum pressure. The actual drop of temperature to maximum pressure 
can be shown to be always less than R8 /2 (&+/), which in this case is 240°, 
a result which justifies the procedure by which, in the calculation of the 
maximum pressure, h has been taken to be a linear function of the 
temperature. 

In this case also the integration of (11) can be carried out analytically, 
6 being given in terms of h from (la), and the value of log(Ui/Uo) can be 
expressed in terms of B and z 0i Ui being the value of U at the instant when 
the charge is all burnt. For practical purposes it would be necessary to 
tabulate log(Ui/U ) as a double-entry table in B and % but various methods 
of approximate integration of (11) can be devised, which are sufficiently 
accurate for the purpose, and which avoid the heavy arithmetical work 
involved in the calculation of such a table. 

.As an example, the theory has been applied to the case of an 11-inch gun, 
using M.D. cordite, for which we may take m\ = 238 lb., M — 760 lb., 
n = 0*15 in. The value of the constant a for the rate of burning is taken 
from Mansell's investigation, and is 0*611 in. 3 /Ton~sec. for the charge at 
60° F. The initial pressures have been taken as 0, 1, and 2 Tons/in. 2 
respectively, and the corresponding indicator diagrams are plotted in fig. 3. 
They show the effect of raising the initial pressure on both the maximum 
pressure and on the position of the projectile when the charge is all burnt. 
It has also, of course, a considerable influence on the total therm o-dynamical 
work done, which is proportional to the area under the curves, its 
magnitude for the three cases being 46,590, 47,710, and 48,380 foot-Tons 
respectively. 

The value of the initial pressure, p , which must be taken in order to maka 
the experimental and theoretical maximum pressures agree, depends naturally 
on the value of the constant a chosen as the standard, and is usually of the 
order of 2 to 4 Tons/ in. 2 . The value of p having been chosen in this way, it 
is found that in the largest guns the experimental kinetic energy of the 
projectile is from 92 to 95 per cent, of the total energy measured in the 
indicator diagram. 

The above statements cannot, of course, apply to all guns and to all 
charges ; they are only approximate, but show the order of magnitude of the 
empirical factors, which in any form of theory have to be introduced in order 
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Fig. 3. 

to take into account the resistance of the driving band and the losses of 
energy in heat and friction. 



vol. c, — A. 



2 i 
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Equations for Zero Initial Pressure. 

The equations of the theory are considerably simplified if it is assumed 
that the projectile begins to move directly any gas has been given off from 
the burning charge. The corresponding equations can easily be found by 
making y zero in (9), (11), (12), and (15), the whole course of the explosion 
then depending on the single parameter B. 

Applied to any particular case, these equations will always give too low a 
maximum pressure, but they may be used with a modified interpretation so 
as to be of practical value. Apart from the constants in the expression for h 
as a function of 6, there are two other constants, # and a, which appear in 
the equations. If it be assumed that the frictional force in the rifling can 
be allowed for by a virtual increase in the mass of the projectile, then it is a 
simple extention to consider all such forces which give rise to losses of energy 
as allowed for in this way, which is equivalent to a decrease in the funda- 
mental constant B. 

Taking the constants <9 &nd a, as determined from theoretical or experi- 
mental considerations, we can allow for all neglected losses of energy by 
finding in what ratio B must be diminished in order, for example, to give the 
correct maximum pressure. The same ratio — a coefficient of ignorance — 
may then be used in calculations for a similar type of gun. 

This method is essentially that employed by Ingalls in his 'Interior 
Ballistics/ Ingalls has two disposable constants, corresponding in some ways 
to the #o and a of this theory, and he determines them separately, so as to 
make the experimental and calculated values of the maximum pressure and 
muzzle velocity agree. These values can then be employed in future calcula- 
tions with the same propellant and a similar type of gun. 



